Superconducting nonequilibrium transport through a weakly interacting quantum dot 



00 
O 

o 

(N 



(N 



S3 
O 

o 

I 

5h 

Oh! 

33 

o3 



c 

o 

> 

o 
in 



o 

oo 
o 



X 



L. Dell' Anna, A. Zazunov, and R. Egger 
Institut fur Theoretische Physik, Heinrich-Heine-Universitat, D-40225 Diisseldorf, Germany 

(Dated: March 26, 2008) 

We study the out-of-equilibrium current through an interacting quantum dot modelled as an 
Anderson impurity contacted by two BCS superconductors held at fixed voltage bias. In order to 
account for multiple Andreev reflections, we develop a Keldysh Green's function scheme perturbative 
in the dot's interaction strength. We find an unexpected enhancement of the current due to repulsive 
interactions for small to intermediate lead-to-dot couplings. 

PACS numbers: 73.63.-b, 74.45. +c, 74.50.+r 



I. INTRODUCTION 



Superconducting transport through low-dimensional 
nanoscale structures is currently attracting consider- 
able interest. Gate-tunable Josephson currents through 
nanowire-based quantum dots have been reported^ and 
similar setups have been realized using (short) carbon 
nanotubes 2 -^ and metallofullerene molecules^ Nonequi- 
librium transport in such systems contacted by super- 
conducting electrodes has been a particular focus of rc- 

mainly caused by an 



cent experimental effort j^iiSiSiZii 
interesting interplay between interaction effects (on the 
quantum dot) and superconducting correlations (due to 
the electrodes). One remarkable consequence is the ob- 
servation of an 'even-odd' effect (as a function of the dot's 
occupation number) in the conductance^^ This effect is 
presumably caused by the absence or presence of Kondo 
correlations. In this paper, we analyze superconducting 
transport through an interacting quantum dot for the 
simplest case of a single spin-degenerate level with repul- 
sive on-site interaction energy U > (Anderson impu- 
rity), contacted by two wide s-wavc BCS superconduct- 
ing electrodes with identical gap A. For simplicity, we as- 
sume that both lcad-to-dot couplings (hybridizations) are 
equal, = = T, and consider the two electrodes held 
at potential difference (voltage bias) V. Under a system- 
atic perturbative expansion in the interaction strength U , 
we compute the I — V characteristics, in particular for 
the interesting subgap regime eV < 2A, where multiple 
Andreev reflection (MAR) processes provide the domi- 
nant transport mechanism. A theory of coherent MAR 
has been originally developed for superconducting point 
contacts^ with the essential assumption that charging 
interaction effects inside the contact can be neglected. 
The problem of resonant MAR through a noninteracting 
quantum level has been treated in Refs. [TTIfT^fl3l . 

While the interacting problem in equilibrium has been 
theoretically studied by many authors^ the correspond- 
ing nonequilibrium problem is more difficult and far 
less understood. Previous approaches can be broadly 
grouped in three classes, (i) By ignoring MAR pro- 
cesses in the Coulomb blockade regime, additional side- 
peaks in the differential conductance at eV — 2A and 
2(A + U) were predicted^ reflecting the singularity of 



the BCS spectral density of the leads, (ii) Different 
mean-field schemes have been proposed ) 12 i 16 ' 17 based on 
slave-boson or Hubbard-Stratonovich-path-integral ap- 
proaches. These calculations predict an overall suppres- 
sion of the current by the interactions. This suppression 
is obtained only for sufficiently repulsive interactions, 
while there is no interaction effect for weak interaction^ 
(iii) A Fermi liquid approach valid in the deep Kondo 
limit has been proposed^ Here we do not discuss the 
Kondo regime, but instead focus on the limit of weak 
interactions, U/T < 1, where a controlled perturba- 
tive expansion in the small parameter U /T is possible. 
Note that this approach still allows for arbitrary ratio 
r/A. Such calculations have been carried out for normal- 
conducting (A = 0) electrodes recently ) 19 : 20 ! 21 and we 
here generalize them to superconducting electrodes. The 
case U < r is of experimental relevance for the under- 
standing of superconducting transport through quantum 
dots or molecules with good lead-to-dot couplings. 

The structure of the remainder of this article is as 
follows. In Sec. [Til we discuss our perturbation theory 
approach to superconducting transport through an An- 
derson dot, and its numerical implementation. Results 
for the current- voltage characteristics are shown and dis- 
cussed in Sec. Mil The appendix contains qualitative 
arguments for the current enhancement found at T < A, 
based on an evaluation of the Josephson current. We 
often set h = e = 1 . 



II. PERTURBATIVE APPROACH TO 
SUPERCONDUCTING TRANSPORT 

We consider the canonical Anderson impurity model, 
H = Hd + Ht + Ht, + Hr, where a single- level dot 
with spinful fermion d a (Ho) is tunnel-coupled (Ht) to 
left/right superconducting reservoirs Hl/r held at chem- 
ical potential difference eV. The isolated dot corresponds 
to (n CT = 4d„ = 0, 1) 

H D = E (n 1 +n l ) + Un 1 n l = e (n T + n l) - 7T ( n T ~ n l) 2 - 

(1) 

The 'noninteracting' model below is taken to contain the 
interaction level shift eo = Eq + U/2 of the bare level 
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Eq. The leads arc described by a pair of s-wave BCS 
Hamiltonians in the standard wide-band limit. We are 
interested in the V ^ case, and take the same real- 
valued gap parameter A > for both electrodes. Using 
the Nambu vector ^J k = (V'j'.fe.T > V 1 ] -k f° r electrons 
in lead j = L/R, we thus have (we put h = e = 1 in 
intermediate steps) 

H J = E *5* (( fc V2m - e F )a, + A<7 a ) tt jfc) (2) 



with Pauli matrices cr,; (t;) in Nambu (Keldysh) space. 
Using the Nambu vector d = (dp d{) T and T = 7rz/o|io| 2 
for (normal) lead density of states vq, the lead-dot cou- 
pling is 

H T = h & jk <Tze ±i °* vt / 2 d + h.c, (3) 

k,j=L/R=± 

where the voltage V enters via the time-dependent phase. 
We now define the Keldysh-Nambu Green's function for 
the dot fcrmion, 

G^,i/)=-i(f c [d a (t s )dUt s! )}}, (4) 

where a, a' = 1,2 (s,s' — 1,2) are Nambu (Keldysh) 
indices, and Tq is the time-ordering operator along the 
Keldysh contour. Accordingly, t s denotes a time taken 
on branch s of the Keldysh contour. It is convenient to 
use the Fourier decomposition 2 ^ 

oo j- J 

G ab (t,t')= Yl y^e-^'+^'G^M, (5) 
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where a,b = 1, . . . , 4 denotes Nambu-Keldysh indices de- 
fined by a = a + 2(s — 1), and ui n = u) + nV for uj within 
the 'fundamental' domain F = [-V/2,V/2\. For fixed 
u> E F, the Dyson equation for the full Green's func- 
tion G (the check notation refers to the Keldysh-Nambu 
structure), 



G~ 



Gn 



E, 



(6) 



then becomes a matrix equation suitable for numerical 
inversion. Here, interaction effects are encoded in the 
sclf-cncrgy E. After integrating out the lead fermion de- 
grees of freedom, the noninteracting Green's function Go 
is 

G (7.nmM = K - ^z)T z 8 nm - V Y Jj,nm(u). (7) 

j=L/R 

The self-energy due to tracing out the respective lead is 
given by the Nambu matrix 

*fj=L/R=±,nm(v) = ( 8 ) 

X{oj n T V/2) 6 Y((j n T V/2) \ 
\ ^m,n±l Y(w n ±V/2) 5 nm X{u n ±V/2) J 



with Keldysh matrices Y(ui) = —AX(uj)/lli and 

f-72fe T - 1^1 <A 

X(u>) = I iH (2f u -l -2f u \ . , 

where f u = 1/(1 + e u) / kBT ) is the Fermi function. The 
steady-state dc current through the left/right junction 
then follows as 



2~tr [o"z7z,/fi,nm( w )G TO nH] 

(9) 

where the trace is only over Nambu space, and (H — ) 
refers to the (12) Keldysh component. Current conser- 
vation, II = Ir = I, is fulfilled for all results below. 



First order 





FIG. 1: Electron-electron interaction self-energy diagrams 
taken into account in this paper: (a) First order (left), (b) 
Second order (right). 

Since the exact self-energy £ is not known, we proceed 
in a perturbative fashion, starting with the first-order 
self-energy in Fig.QJa), made self-consistent by using the 
full G in the diagram. It is convenient to introduce the 
four-point vertex, cp. also Ref. 0, 



A„ 



bed 



2 (SabTcd + ScdTab ~ S a dT c b ~ ScbTad) , (10) 



where r = ctqt z = diag(l, 1, —1, —1). Only 8 out of the 
possible 256 entries of the tensor A are nonzero, with 
values ±1. With this convention, the complete first-order 
self-energy is given by 



-,(1)06 



iUA a 



bed 



"2T 



G n ' t n'+ m {^ )j 



(11) 



where we use the sum convention for c, d. Note that 
this self-energy is independent of (n, u>), but still depends 
on the 'off-diagonal' frequency index m. In time repre- 
sentation, to 7^ contributions come with phase factors 
e ±2imVt anc j fang correspond to anomalous (pairing) cor- 
relations. The presence of the off-diagonal harmonics in 
the self-energy is a consequence of the coherent MAR 
transport regime considered here. The role of coherence 
is particularly important for the interplay between MAR 
processes and charging effects in a quantum dot with a 
relatively strong coupling to the leads. 

Let us at that stage briefly compare the first-order self- 
consistent approach based on Eqs. ^ and (jTTJ) to the 
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mean-field approximation of Ref. [lj| The latter effec- 
tively considers only the time averaged components of 
the self-energy (fTTj) corresponding to the two first terms 
of Eq. (jTUJ) , thereby discarding all harmonics with m 7^ 
and exchange terms in Eq. . The resulting self-energy 
contributions taken into account in Ref. correspond to 
7+ To and 7_r z . The scalar constants j± can be written 
as 

7+ = ^Ylf ^( T ^nv(w')), 

n' ^ 

7- = 4E/ ^Tr(G n , n ,(u/)), (12) 

where the trace is over both Nambu and Keldysh spaces. 
Under this approximation, there is no interaction effect 
on the current I(V) below some critical value U c . In fact, 
nontrivial stable solutions 7± 7^ for the self-consistency 
equation (fT2|) exist only for U > U c ^ where U c depends 
on V, r and A. The symmetry-broken phase with 7_ 7^ 
corresponds to a spin-polarized dot, and one then finds a 
Coulomb blockade suppression of the current^ However, 
for normal leads, serious problems with spin-polarized 
out-of-equilibrium mean-field solutions for the Anderson 
dot have been identified recently^ and those arguments 
also apply to the superconducting case. A typical value 
is U c « T for V « T as A/2, and we shall always limit 
ourselves to U < U c where no such problems arise. In our 
calculations, the actual value of U c follows from the nu- 
merical solution of the self-consistency problem, and we 
can thereby ensure that no spin-polarized solutions are 
present. In contrast to the mean-field scheme of Ref. [ill 
however, the full first-order self-consistent approach gen- 
erates sizeable interaction corrections even for small U, 
see Sec. IIIIl These corrections are not just a matter of 
numerical accuracy but reflect the importance of on-dot 
pairing terms. 



B. Second order 

To go beyond the self-consistent first-order approxima- 
tion given by Eq. (fTTj) , we have to evaluate the second- 
order diagram shown in Fig. [TJb). -Due to the large nu- 
merical effort in evaluating this diagram, we restrict our- 
selves to a non-selfconsistent scheme at this point, i.e., we 
use the solution G of the first-order problem to evaluate 
Y,( 2 \ As is well known ; 24 ' 25 under such a scheme current 
conservation is only ensured for the particle-hole sym- 
metric case, eo = 0, see also Ref. [22|. We therefore show 
second-order results only for e = 0. The second-order 
Nambu-Keldysh self-energy reads in time representation 
(sum convention) 

£ (2)a6 (M') = ^KfghKbcd (13) 
x G fe (t,t')G d9 {t' \t)G hc (t,t'). 



To avoid numerically expensive frequency convolutions, 
it is convenient to first compute the Green's function in 
time representation according to Eq. (|5|), then evaluate 
the self-energy in Eq. (fT5|) . and finally transform this 
result back to frequency space to use it in the Dyson 
equation ©. Notice that Eq. (fT3"|) corresponding to the 
skeleton diagram in Fig. [TJb) represents only a part of 
all possible second-order contributions. The rest is given 
by the time-local piece 

~U 2 A abcd A efgh J dt" G de (t,t")G h v(t",t")Gf c (t",t), 

which has already been taken into account by our self- 
consistent first-order solution. 



C. Calculation of current 

The numerical implementation of the above perturba- 
tive approach is straightforward. To evaluate the current 
I(V) from Eq. ©, we partition the frequency summa- 
tions into windows of width V, and impose a bandwidth 
cutoff uj c , such that \ui n \ < lu c . In our calculations, we 
use uj c = 10A, but the precise choice is not critical. 
We then discretize the fundamental frequency domain 
F with a step-size 8cu, and use a fast Fourier transform 
routine to switch between time and frequency represen- 
tations. (The efficient evaluation of the second-order self- 
energy requires to employ the time representation, while 
the Dyson equation needs the frequency representation.) 
Typically, we found Suj = 0.005A sufficient for conver- 
gence. The matrix inversion in Eq. ([6]) is then performed 
for each ui £ F separately, involving matrix dimensions 
of the order of o; c /|V|. We refer to Ref. I22I for further de- 
tails of the numerical implementation in the related case 
of a phonon-mediated interaction. 

In a first step, we solve the first-order self-consistent 
problem posed by Eqs. ([6|), ([7j) and (fTT|) . This solu- 
tion proceeds iteratively, where the stability or instabil- 
ity of the solution for G is checked carefully by probing 
small deviations around it. The iterative solution can 
in fact be carried out with very modest computational 
effort, and quickly converges to a unique solution (as 
long as U < U c ). In a second step, we then use this 
converged first-order Green's function to evaluate S*- 2 -* 
according to Eq. () 13ft . and to finally compute the cur- 
rent from Eq. ([9]). The second-order calculation is quite 
time-consuming for low bias voltage, where many MAR 
orders need to be taken into account, and we therefore re- 
stricted our calculations to eV/A > 0.2. As consistency 
check for our numerical code, we have reproduced known 
results for U = 0, see Refs. [lllfl2lf l3l and the correspond- 
ing perturbative-in-f/ results for normal-conducting leads 
(A = 0), see Refs. H9II20I W e have also reproduced the 
respective results of Ref. [ly when implementing their ap- 
proximations. As additional check, Green's function sum 
rules, such as tr [r z a z G(t, t)] = at coinciding times, 
have been verified. 
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III. RESULTS AND DISCUSSION 

Next we discuss numerical results obtained under the 
perturbative approach described in Sec. [TTJ All results 
are for T = 0, and unless noted otherwise, we put 
U/T = 0.5 which is sufficiently small to ensure U < U c 
for all investigated V/A and T/A but large enough to 
produce significant interaction corrections to the I — V 
characteristics. We will focus on the most interesting 
subgap regime. eV < 2A. The excess current I exc = 
limy^ 00 [/(y, A) — I(V, A = 0)] has also been computed. 
The interaction contribution dl exc to this quantity turns 
out to be generally small, similar to what is found for 
the case of phonon-mcdiatcd interactions^ Remarkably, 
this interaction correction is positive for T < A, point- 
ing towards a current enhancement. This trend is quite 
generic and discussed next. 



0.15 



0.06 ; 




FIG. 2: Interaction correction to the current (currents are 
always plotted in units of e A / h) from the self-consistent first- 
order approach, for U/T = T/A = 0.5 and various eo/A. The 
inverse voltage scale is taken to compare with standard MAR 
features. Inset: Full I — V curves for same parameters. 

Let us start by showing results obtained from the first- 
order self-consistent scheme (i.e., without the second- 
order self-energy). In that case, by virtue of self- 
consistency, we have the freedom to vary eo without spoil- 
ing current conservation! 24 i 25 Representative numerical 
results for the voltage-dependent interaction correction 
to the current, I(U)-I(U = 0), are shown for T/A = 0.5 
in Fig. [2] For all results shown here, we have U < U c , and 
the approach of Ref. [lfl would not yield any interaction 
correction. However, we find significant interaction ef- 
fects for U < U c within the full first-order self-consistent 
approach. These effects arc due to the time-dependent 
(in 7^ 0) parts of the self-energy (fTTj) . which contain pair- 
ing order parameters on the dot. For instance, at the 
symmetric point eo ~ 0, we find by perturbation theory 
in U that j± = 0, but the complex- valued m = 1 pairing 
term 




-0.02 



FIG. 3: Same as Fig. [2] but for eo = 0. The dashed curve 
gives the first-order self-consistent result, while the solid curve 
includes also the second-order contribution. 



r=? 




.(1),12 



(14) 



FIG. 4: Same as Fig. E] but for T = 2A. Inset: I — V curve 
from second-order perturbation theory and for U = 0. 



stays finite. This ^-independent off-diagonal Nambu 
component of the self-energy, absent in the normal 
(A = 0) case, describes the effect of interactions on the 
proximity-induced pairing correlation on the dot. At the 
mean-field level, 7_ is dominant for large U, while terms 
like S dominate for small U. Similar contributions with 
|m| > 1 exist and are kept in our self-consistent first- 
order calculations, but they turn out to be significantly 
smaller. 

Quite remarkably we find I(U) > I(U = 0) for most 
voltages and/or dot level energies eo, pointing to an en- 
hancement of the MAR-mcdiated current by repulsive in- 
teractions. We have persistently found this unexpected 
feature throughout the parameter regime T < A, also 
when including the second-order self-energy, see below. 
A similar (but weaker) enhancement can be found ana- 
lytically for the critical Josephson current of this system 
in equilibrium, see Appendix. The current enhancement 
is reminiscent yet different from the 'antiblockade' behav- 
ior due to dynamical Coulomb blockade effects on MAR 
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transport discussed in Ref. [26|. It is also consistent with 
the crossover from current enhancement to decrease with 
growing T for phonon-mcdiatcd interactions and normal- 
conducting leads ^ 

To illustrate the role of the second-order contribution 
for U/T = 0.5, we now focus on the symmetric case 
eg = 0, first taking again T/A — 0.5. The results of 
the first- and second-order calculations are compared in 
Fig. [3l Notice that the second-order correction, which 
is the leading time-nonlocal term in the perturbative ex- 
pansion, becomes more and more important when lower- 
ing the voltage. In agreement with the conclusion drawn 
from the first-order self-consistent calculation shown in 
Fig. [2j a clear enhancement of the current by interactions 
can be observed for a broad range of voltages. This en- 
hancement is especially pronounced for voltages slightly 
below the odd MAR peaks located at eV = 2A/(2n+ lj. 
As indicated by our results for larger T, see Fig. 0] for the 
case T/A = 2, the interaction-induced enhancement of 
the current is restricted to small T/A. For larger T/A, 
the current instead is weakly suppressed by interactions. 
The same pronounced MAR peak structure as in Fig. [3] 
can be observed in the pair order parameter S defined in 
Eq. (fT4j) . whose absolute value is shown in Fig. [5] The 
fact that the characteristic MAR features still appear at 
eV = 2A/n for the interacting case indicates that, at 
least for small U, the number of Andreev reflections is 
not affected by the Coulomb interaction. This is in con- 
trast to the inelastic MAR picture for phonon-mediatcd 
interactions, as discussed in Ref. [U Moreover, our re- 
sults indicate that the Andreev quasircsonance a 1 1 ' 13 are 
not shifted away from the gap subharmonics. 
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FIG. 5: Absolute value of 5, see Eq. ([14]) . in units of A, for 
the parameters in Fig. [3] 

As follows from our numerical analysis, the magnitude 
of the difference current is mainly determined by the first 
harmonics of interaction- mediated pairing, Eq. (|14[) . and 
can be roughly approximated as 

I(U)-I(Q)~(e/H)\8\A/T. (15) 

This can be seen, for instance, from comparison of the 



curves in Figs. [5] and [2] (for eo = 0). A similar expres- 
sion for the difference current, Eq. (fT5"]) . is also obtained 
from a simple Fermi golden rule calculation, by analyz- 
ing the dynamics of Andreev (subgap) states at very low 
voltages, eV -C A. The corresponding correction to the 
transition rate from Andreev states into the continuum is 
determined by the imaginary part of the Andreev state 
sclf-encrgy ^^(w) = 5\J A 2 — lo 2 /T, while the total es- 
cape probability leading to the difference current is given 
by an expression similar to Eq. (3) in Ref. [l8|. 

The above results also suggest that as a function of 
the ratio T/A, there should be a crossover from enhance- 
ment to suppression of the current around T/A ps 1. 
This is what we get when fixing the voltage and chang- 
ing T/A. In Fig.0 we have chosen V = 0.6A, where the 
self-consistent first-order approximation gives the main 
contribution, and plotted I(U) — 1(0) either for fixed 
U/T = 0.5 (solid line), or for fixed U/A = 0.25 (dashed 
line). The two curves both cross zero approximately at 
the same value, T ~ A. 




T/A 

FIG. 6: I(U) - 1(0) versus T/A at V = 0.6A, e = 0, fixing 
U/T = 0.5 (solid line) or U/A = 0.25 (dashed line), from 
self-consistent first-order perturbation theory. 

In conclusion, we have presented a theory exploring 
the effect of weak interactions on superconducting trans- 
port through a quantum dot. Employing second-order 
perturbation theory, valid for U < T, we find an un- 
expected enhancement of the subgap current against its 
noninteracting value when the hybridization T is smaller 
than the BCS gap parameter A. The perturbation the- 
ory scheme pursued in this paper offers controlled results 
in one corner of the parameter regime, and in contrast to 
previous mean-field theories, we predict significant inter- 
action corrections even for weak interactions. 
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APPENDIX A: JOSEPHSON CURRENT 

In this appendix, we briefly show that the interaction- 
induced current enhancement found in the I — V curves 
for T < A discussed in Sec. HIT] also appears in the equi- 
librium Josephson current-phase relation (where <j) is the 
phase difference across the dot) for the same model. We 
consider the corresponding first-order self-consistent the- 
ory in equilibrium, for simplicity at eo = only. For 
small U, no polarization is present, 7 = 0, and only 
the proximity- induced mean- field parameter 5 = U (cijdj) 
gives an effect. Assuming real- valued S, the T = self- 
consistency equation reads 



U 



2tt 



+ - sy 



(Al) 



where 



= U) 1 



A 5 



/3a, = 



TAcosO/2) 
V^ 2 + A 2 



The Josephson current is then given as 



erA 
1= —- sin ^/2 
irn 



0u 



, Vw 2 + A 2 al + - 5)' 

(A2) 

The presence of S in Eq. (|A2|) generally causes two coun- 
teracting effects: there is a decrease of / due to the nu- 
merator, but an increase due to the appearance of 8 in 
the denominator. Which of these is more important can 
only be clarified by detailed calculation. We present ana- 
lytical evaluations valid for \5\ <C T| cos(</>/2)|, separately 
for the regimes T/A < 1 and F/A > 1. 

Let us first discuss T/A <C 1, where Eq. (|A1[) yields 



u 



2|rcos(>/2) - S\ 
2rcos(>/2)cos 



(A3) 



_1 | rcos(0/2)-<5 I 
A I 



A 2 - [rcos(^/2) - 8} 



The interaction correction to the Josephson current now 
follows from Eq. ljA2|) . 



Sl~ ^tan(4>/2) Q/x^sgncos^)-^ 



with 



l + 2(r/A)|cos(<ft/2)| 
[1 + (T/A)|cos(0/2)|] 2 



= i + o(r 2 /A 2 ). 



In the extreme limit T/A -» 0, Eq. (|A3j) yields 5 = 
(C//2)sgncos(0/2), and then 5/ = 0. Inspection of 
Eq. (fAltj) for finite T/A <C 1 shows however that \S\ < 
U /2 for <f> ^ tt. As a result, the interaction current to the 
Josephson current for T <C A turns out to be positive, 
albeit numerically small. Using Eq. (|A3[) we find 

eU r 

51 ~ — — — sin(0/2) sgncos (0/2). 

We believe that this effect is related to the enhancement 
of the current at finite bias V in the regime T <C A 
discussed in Sec. Mil 

On the other hand, for T > A, Eq. ([AT]) is solved 
by S ~ E/Aln(r/A)cos(0/2)/(7rr), and the lowest-order 
interaction correction to the Josephson current is 



SI ~ ^tan(0/2) ("A /2( ^ sgncos(0/ 2)- A 



with f 2 {4>) = 1 + cos 2 (0/2)/4. Due to the large ln(r/A) 
factor appearing now in 8, the Josephson current will in 
general be decreased by interactions for r> A. 

We therefore find the same qualitative picture as for 
the nonequilibrium current in Sec. IIIII The equilibrium 
Josephson current can also be slightly increased by weak 
repulsive interactions for weak hybridization, T/A <C 1, 
but is decreased in the opposite limit. However, the in- 
crease of the Josephson current for weak hybridization 
turns out to be much smaller than for the corresponding 
nonequilibrium current. 
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